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Abstract. We have identified 22 new variable red giants in 47 Tuc and determined periods for another 8 previously known 
variables. All red giants redder than V-/c = 1.8 are variable at the limits of our detection threshold, which corresponds to 
SV ~0.1 mag. This colour limit corresponds to a luminosity log L/Lq=3. 15 and it is considerably below the tip of the RGB at 
logL/Lo=3.35. Linear non-adiabatic models without mass loss on the giant branch can not reproduce the observed PL laws 
for the low amplitude pulsators. Models that have undergone mass loss do reproduce the observed PL relations and they show 
that mass loss of the order of 0.3 M© occurs along the RGB and AGB. The linear pulsation periods do not agree well with the 
observed periods of the large amplitude Mira variables, which pulsate in the fundamental mode. The solution to this problem 
appears to be that the nonlinear pulsation periods in these low mass stars are considerably shorter than the linear pulsation 
periods due to a rearrangement of stellar structure caused by the pulsation. Both observations and theory show that stars 
evolve up the RGB and first part of the AGB pulsating in low order overtone modes, then switch to fundamental mode at high 
luminosities. 

Key words, stars: late-type - stars: AGB and post-AGB - stars: oscillations - stars: mass-loss 



1. Introduction 

It is now well established that pulsating red giant stars lie on 
a series of up to six parallel period-luminosity (PL) sequences 
(e.g. Wood et al. lT99^ Wood l2000l Ita et al. l2n^ Soszyiiski et 
al. 120041 Fraser et al. 12005 > . The stars known to populate these 
sequences are generally field variables so that their metallici- 
ties, masses and ages are not known individually, although the 
luminosities are well known because the stars lie in stellar sys- 
tems at a known distance, such as the SMC and LMC. The dis- 
persion in the PL relations is much larger than the amplitude of 
the pulsation, especially for the small amplitude variables, and 
this dispersion is almost certainly a consequence of the disper- 
sion in mass and metallicity at a given pulsation period. 

Comparison of the PL relations with theory (Wood et al. 
I1999> shows that four of the PL sequences (sequences A, B 
and C of Wood et al. 119991 or sequences A, B, C and C of 
Ita et al. l2004> can be broadly explained by radial pulsation in 
the lower order modes. However, a real test of the theoretical 
models requires comparison with stars of known metallicity, 
age and (initial) mass since the fits of theory to data in Wood et 
al. (I1999> and Ita et al. (I2004> are not particularly good. 

Globular clusters are well suited for carrying out such 
a comparison, and for studying the relation between pulsa- 
tion, mass loss and stellar evolution along the giant branch. 



since fundamental parameters like initial mass, luminosity and 
metallicity are well known. However, only a small number of 
pulsating red giants, hereinafter referred to as long-period vari- 
ables or LPVs, is known in any given globular cluster (Clement 
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et al. l2001t . with 47 Tuc having the most known LPVs (14). 
This is an insufficient number of variables for clearly defining 
PL sequences, especially as many of the periods of the LPVs 
are poorly determined. The origin of this problem is that most 
search programs for variable stars in globular clusters were de- 
signed to optimize detection and period determination for vari- 
ables with periods of about one day or less. LPVs with their 
periods of a few ten to a few hundred days have therefore not 
been well surveyed and the existing samples are far from com- 
plete. 

To rectify this situation, we have started a search pro- 
gram for long period variables in Galactic Globular clusters 
(Lebzelter et al. .2004j . The relatively large number of already 
known or proposed long period variables in 47 Tuc made it a 
good starting point. Here, we present the results for 47 Tuc 
(NGC 104), the first cluster analyzed. We adopt the follow- 
ing properties for 47 Tuc: {m - M)y=13.5+0.08 (Gratton et al. 
EnnS; metallicity [Fe/H]=-0.66 (Carretta & Gratton llWTt : in- 
terstellar reddening E(B - y)=0.024; an age of 11.2+1.1 Gyr 
(Gratton et al. l2003> : and a turnoff mass between 0.86 and 
0.9 Mq. In a recent paper discussing radial velocity variations 
in known LPVs in 47 Tuc (Lebzelter et al. l2005t we gave more 
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discussion of the properties of 47 Tuc and its variables and we 
refer the reader to that paper for details. 

2. Observations 

2.1. Mount Stromlo data 

Our monitoring program of 47 Tuc started in August 2002 
at the 50 inch telescope at Mount Stromlo. This telescope 
was equipped with a two channel camera used earlier for the 
MACHO experiment (Alcock et al. ll992> . The camera ob- 
tained two images in two broad band ranges at the same time 
(Marshall II 994> . These passbands did not correspond to stan- 
dard filters but the blue one had a mean wavelength similar to 
Johnson V. The camera covers a field of about 0.5 square de- 
gree on the sky with a pixel scale of 0.62 arcseconds. The whole 
cluster could thus be observed with one observation, and most 
of the cluster was on one CCD of the MACHO camera (there 
are four CCDs per filter band). Observations were obtained 
once to twice a week. A few frames were lost due to techni- 
cal problems or bad seeing conditions. Our monitoring came 
to an early end after about five months when Mount Stromlo 
Observatory was destroyed by a bush fire. Altogether we col- 
lected 15 useable frames over this time span. All observations 
were done in queue observing mode. 

For the determination of light curves and the detection of 
variables, only the blue frames were used. The light amplitude 
of long period variables is typically larger in the blue than in the 
red (e.g. Fox ll982t making the detection of variables and the 
determination of periods easier Furthermore, many bright stars 
on the red frame were over-exposed and thus not useable. The 
blue detector had an area of bad pixels as well as a few scat- 
tered dead pixels. Consequently, due to small positional shifts 
between the different observing nights some stars could not be 
measured on all 15 frames. 

2.2. CTIO data 

From August 2003 to January 2004, we continued the monitor- 
ing with ANDICAM at CTIO's 1 .3m telescope operated by the 
SMARTS consortium. The size of the CCD field is 6x6 arcmin- 
utes and the pixel scale is 0.37 arcsec/pixel. A description of the 
camera is given in DePoy et al. ( t2003i) . Due to the smaller FOV 
we had to make a mosaic of images to cover most of the clus- 
ter Even then, we were limited to the central part, and some of 
the outer areas covered by the Mount Stromlo data could not 
be observed. Observations were done in V and /c. Observations 
were scheduled roughly once per week giving a total of 16 use- 
able frames. As for the Mount Stromlo data, observations were 
done in queue observing mode. 

2.3. Other data 

Six observations of the cluster in V and 1^ were taken in ser- 
vice mode with the WFI at the ESO 2.2m. These observations 
were obtained between June and September 2002. Four addi- 
tional observations taken in July/August 2003 were kindly pro- 
vided by Laszlo Kiss and collaborators. These observations in 



V and /c were obtained with the Siding Spring Observatory 40- 
inch telescope, helping to reduce the gap between the Mount 
Stromlo and the CTIO data set. 

For a different project (Lebzelter et al. l2005> . we did a 
short photometric monitoring of some parts of 47 Tuc with 
ANDICAM at CTIO's 1.3m telescope in March to May 
2002. These data were obtained in the V-filter with the old 
ANDICAM CCD which had severe Hmitations (only half of 
the chip was useable). These data were used only in rare cases. 

3. Data reduction 

For the Mount Stromlo and the CTIO data, flatfield and bias 
correction was done as part of the standard data pipeline. The 
40inch data were reduced applying standard data reduction 
with MIDAS, while the WFI data were reduced with the cor- 
responding IRAF package. For the detection of variables and 
the measurement of the light curves on the Mount Stromlo 
and CTIO frames we used the image subtraction code ISIS 2. 1 
by Alard J2000t . First, the two data sets were analyzed sep- 
arately. The reference flux of the identified variables required 
to produce light curves from the image subtraction measure- 
ments was derived using the PSF fitting software written by 
Ch. Alard for the DENIS project. A description of this code 
can be found in Schuller et al. ( I2003t . For the photometric cali- 
bration of the CTIO data we used standard stars from Landolt's 
field RU 149 (Landolt ll992l l. Photometric accuracy of the re- 
sulting light curves could be analyzed by comparing measure- 
ments of variables in overlapping parts of the CTIO mosaic. 
Typical deviations were of the order of 0.007 mag. We then 
used about 40 non-variable cluster stars ranging in V - / be- 
tween 0.9 and 2.1 to link the Mount Stromlo measurements to 
those from CTIO. Typical errors in the Mount Stromlo magni- 
tudes were of the order of 0.013 mag. 

For various reasons (field of view, pixel scale, damaged 
parts of the CCD, depth of the observation) not all variables 
could be measured on both the CTIO and the Mount Stromlo 
data. In these cases only one of the data sets could be used for 
the analysis. 

Our aim was to detect and measure long period variables 
on the upper part of the giant branch only. Thus we selected for 
analysis only those variables that varied on time scales of more 
than approximately 30 days with a total light amplitude of at 
least 0.1 mag in V (or the blue MACHO pass band). For stars 
with shorter periods, our sampling was not sufficient to esti- 
mate a useful period. We searched for periods in the selected 
light curves using the Fourier analysis code Period98 (Sperl 
II998t . A maximum of two periods was derived for each star. 
First, the two data sets from Mount Stromlo and CTIO were an- 
alyzed separately. The observations from ESO and SSO, mea- 
sured using PSF fitting, were inserted in the combined light 
curves. Then a period search was done on the whole light curve. 
Generally, the agreement between the periods derived from 
these three data sets agreed very well. Naturally, stars with pe- 
riods exceeding the length of one of the data sets or with high 
irregularity in their light curve lead to deviating results. The 
uncertainties of the periods calculated from a least square fit to 
the complete data sets are in most cases less than a few per- 
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Fig. 1. Light curves of the variables of our survey, except those stars akeady discussed in Lebzelter et al. iZQQSi . Fourier fits with 
up to two periods are shown. Note that both axis scales vary from star to star For A19, LW5, LW7, LW13 and LW14, a break 
has been introduced into the time axis for a better presentation of the light change. Errors on data points before JD2452700 are 
estimated to be 0.013 mag. (Mount Stromlo data) while after this date the errors are estimated to be 0.007 mag. (CTIO data). 
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cent. Of course the semiregular nature of most of the stars in 
our sample means that the derived period is representative only 
for the current variability behaviour. For the long period cases, 
we used only the period derived from the combined light curve. 
Some stars were classified as irregular if no consistent period 
could be found. These stars will not be discussed further in this 
paper 

4. The variables 

Periods have previously been published in the literature for the 
LPVs V1-V8, Vll, V1 3, V1 8 and V21 (we use the nomen- 
clature of Clement et al. l2001> . Light curves for these stars de- 
rived from the current observations are given in Lebzelter et 
al. J2005> . along with some discussion of the origin and re- 
liability of the published periods. These periods are given in 
column 7 of Table[2 Beside these stars, a few more red gi- 
ants have been identified as variables before, namely V15, V16, 
V17, V19, V20, V22, V23, V25, V27, V28 and A19. For some 
of these (V15 to VI 7, V25 and V28) periods were given by 
Fox ( I1982> . A19 was reported variable by Lloyd-Evans (I1974> . 
but no period was given. Later variability surveys, like those 
of Kaluzny et al. ( fT998b and Weldrake et al. ( l2^4b . have fo- 
cused on binaries and RR Lyr stars and thus typically were not 
suitable for determining periods of the stars on the upper gi- 
ant branch due to overexposure. The long period variables they 
report are most likely members of the Small Magellanic Cloud. 

Table[n lists the variables detected or characterized in the 
course of this study. Our sample includes 22 variables detected 
for the first time, and we give the first period determinations 
for six further stars previously known to be variable. All these 
periods are Usted in column 6 of Table[l] We do not include in 
Table^any previously known variables that we were not able to 
monitor (because they fell on CCD defects or outside our field 
of view). As nomenclature for new variables, we use LWxx 
with the numbering going from east to west. Near infrared J 
and K magnitudes were extracted from the 2MASS database 
and they are Usted in columns 4 and 5. The 2MASS coordinates 
are given in columns 2 and 3. 

Light variation versus time are presented in FigurefO to- 
gether with the fits from Fourier analysis. For selected variables 
with a reliable period determination and no long secondary pe- 
riod we show the light curve plotted against phase in Figure|3] 
In Fig. [2 obvious deviations of the data from the fit illustrate 
the semiregular nature of many stars of the sample. It can be 
seen that quite a large fraction of variables shows a long sec- 
ondary period lasting several hundred to several thousand days. 
As these periods exceed the length of the time span monitored, 
we could not determine their length accurately and thus we do 
not give them in Tabled They are included in the fits for illus- 
trative purposes only. The formal ratio between the long and 
the short period in the plot ranges between 4 (LW4) and 100 
(LW15). In most cases the second period is 10 or more times 
longer than the shorter period. 

The Fourier fits to the stars V20, LWIO, LWll, LW12 and 
LW13 deserve special comment. V20 and LWIO both appear 
to have periods of about 220-230 days but, because of the 
time gap in our observations, these periods are quite uncertain. 



Figs.^andl^lsuggest that these two stars have light curves with 
bumps during rising light, a feature found in the light curves of 
several miras in the solar neighbourhood. For LWIO the semi- 
period shows up clearly in the Fourier analysis although this is 
not the case for V20. The fit of LWl 1 is rather bad and the pe- 
riod determined is very uncertain. In LWl 2, the period and am- 
plitude of the light curve clearly changed between the Mount 
Stromlo and CTIO observations and we use separate Fourier 
fits to the two sets of observations. The Mount Stromlo data 
set gave a period of 61 d while the CTIO data give 116 d. The 
fit to LWl 3 is poor, and the plotted fit curves include a phase 
change of 20 days between the Mount Stromlo and CTIO ob- 
servations. Finally, we note that the 2MASS observations could 
not separate LWl 3 from its close companion. 

Multiperiodicity is a well known phenomenon among red 
variables in the solar neighbourhood (e.g. Percy et al. l2003l 
Kiss et al. lI999> . Long secondary periods are well known 
among long period variables, but a physical explanation for 
them has not yet been found (Wood, Olivier & Kawaler [2004t . 

In Figure|3l we show the position of the variables in the 
colour-magnitude diagram using magnitudes from the CTIO 
images or, for a few stars, an alternative source given below. 
Only the giant branches and the horizontal branch are shown. 
Beside the variables listed in Tabled we also marked stars de- 
tected as variables but slightly below our amplitude cutoff crite- 
rion and stars with large irregularities in the light curve so that 
a period could not be determined. All V and V - 1^ values of 
the variables are mean values derived from our light curves. For 
V5, V7, VI 1, V13 and V18 no CTIO data ha ve been obtained. 
Instead, we give the values from Fox J1982> . V3 was also not 
observed at CTIO, but V and Ik light curves were obtained by 
Eggen (II975> and we use the mean magnitudes from his obser- 
vations, converting V - Ik values ioV - Ic using the transform 
in Bessell jl979l l. V2 and LWIO are missing because the stars 
were saturated or unmeasurable on some of our / frames. 

Figure|3]illustrates that all variables in our sample, with the 
possible exception of VI 9, are almost certainly on the upper gi- 
ant branch of 47 Tuc (rather than in the SMC or Galactic halo). 
V19 is located above and to the left of the 47 Tuc AGB. There is 
no indication for a long secondary period in this star (see Fig.[0 
which may lead to a non-representative measurement'. At the 
blue end of the giant branch where V19 is found, its amplitude 
of more than 0.8 mag is rather untypical. Possible explanations 
may be that it has a very close blue neighbour that was not 
separated from V19 in our images, or it may be a star in the 
Galactic halo. 

The most interesting aspect of Figure|3 is that it clearly 
shows that all stars on the giant branch of 47 Tuc redder than 
y - /c ~ 1.8 are variable. 

Finally, we comment on the connection between LPV pul- 
sation and mass loss in 47 Tuc. Origlia et al. J2002> listed five 
stars in 47 Tuc with a high infrared excess, indicating that they 
are surrounded by some dust. All five are located within our 
investigated field. Identification is not simple as Origlia et al. 

' Long secondary periods may be responsible for some of the scat- 
ter along the giant branch since the mean brightness is calculated over 
only part of the light cycle corresponding to these periods. 
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Table 1. The LPVs of 47 Tuc. 



Name 


a (2000) 


6 (2000) 


V 


V-/, 


J 


K 


P[d] 


p 


Remark 


VI 


00 24 12.4 


-72 06 39 


13.15 


3.86 


7.45 


6.21 


221 


in 


1 


V2 


00 24 18.4 


-72 07 59 


_ 


_ 


7.52 


6.29 


203 


203 


1,2 


V3 


00 25 15.9 


-72 03 54 


12.63 


3.54 


7.49 


6.27 


192 


192 


1 


V4 


00 24 00.3 


-72 07 26 


12.34 


2.62 


7.87 


6.69 


165 


82, 165 


1,3 


V5 


00 25 03.7 


-72 09 31 


11.80 


1.99 


8.65 


7.47 


50 


60 


1,4 


V6 


00 24 25.5 


-72 06 30 


11.74 


1.94 


8.54 


7.43 


48 


48 


1 


V7 


00 25 20.6 


-72 06 40 


11.83 


2.30 


8.18 


6.97 


52 


52 


1,4 


V8 


00 24 08.3 


-72 03 54 


12.01 


2.39 


7.94 


6.70 


155 


155 


1 


Vll 


00 25 09.0 


-72 02 17 


12.03 


2.70 


7.91 


6.71 


160: 


52, 100 


4,5 


V13 


00 22 58.3 


-72 06 56 


12.36 


2.23 


8.79 


7.70 


40 + long 


40 


1,4,8 


V16 


00 25 23.2 


-72 1 1 05 


11.65 


2.0 


8.35 


7.23 


41, 88: 


45 




V18 


00 25 09.2 


-72 02 39 


11.67 


1.95 


8.59 


7.47 


83: 




1,4 


V19 


00 24 14.8 


-72 04 44 


11.46 


1.70 


8.79 


7.61 


83 






V20 


00 24 14.5 


-72 05 09 


12.30 


2.53 


8.08 


6.94 


232: 




11 


V21 


00 23 50.1 


-72 05 50 


12.41 


2.85 


8.07 


6.78 


76: + long 




1 


V22 


00 24 08.9 


-72 03 00 


11.80 


1.99 


8.38 


7.18 


62 




8 


V23 


00 24 29.5 


-72 09 08 


11.77 


2.08 


8.68 


7.53 


52 + long 






V25 


00 23 58.9 


-72 02 35 


11.96 


2.50 


8.20 


7.03 


44 


42 


9 


V27 


00 24 15.2 


-72 04 36 


12.11 


2.52 


7.97 


6.77 


69 + long 






A19 


00 24 21.8 


-72 04 13 


12.18 


2.57 


8.08 


6.79 


60 






LWl 


00 23 22.3 


-72 05 40 


11.83 


2.08 


8.39 


7.19 


39 




9 


LW2 


00 23 29.2 


-72 06 20 


11.90 


2.23 


8.26 


7.05 


60 + long 




9 


LW3 


00 23 47.4 


-72 06 53 


11.92 


2.35 


8.13 


6.89 


107 






LW4 


00 23 51.3 


-72 03 49 


11.80 


1.98 


8.56 


7.37 


32 + long: 




9 


LW5 


00 23 53.2 


-72 04 16 


11.74 


1.96 


8.53 


7.32 


40 






LW6 


00 23 54.7 


-72 03 39 


11.85 


2.00 


8.51 


7.30 


41 




9 


LW7 


00 23 56.9 


-72 05 33 


11.66 


1.86 


8.26 


7.18 


42 




10 


LW8 


00 23 57.7 


-72 05 30 


11.82 


1.88 


8.10 


7.12 


27: + long 




9 


LW9 


00 23 58.2 


-72 05 49 


12.34 


2.77 


7.97 


6.74 


74 




9 


LWIO 


00 24 02.6 


-72 05 07 


12.22 


_ 


7.57 


6.40 


110:, 221: 




2,10,11 


LWll 


00 24 03.2 


-72 04 51 


11.95 


1.79 


8.47 


7.40 


36: 




11 


LW12 


00 24 04.0 


-72 05 10 


11.96 


2.35 


8.12 


6.89 


61, 116 




11 


LW13 


00 24 07.9: 


-72 04 32: 


12.56 


2.93 


7.44: 


6.25: 


65: 




6,11 


LW14 


00 24 09.4: 


-72 04 49: 


11.71 


1.99 


8.42 


7.39 


50 




7,10 


LW15 


00 24 11.2 


-72 05 09 


11.76 


2.05 


8.32 


7.13 


46 + long 




9 


LW16 


00 24 13.6 


-72 04 52 


11.88 


1.89 


8.39 


7.32 


29 + long 




9 


LW17 


00 24 16.3 


-72 01 31 


11.81 


1.99 


8.43 


7.20 


41 




9 


LW18 


00 24 20.5 


-72 04 50 


11.91 


2.33 


7.95 


6.77 


65 + long 






LW19 


00 24 23.2 


-72 04 23 


11.85 


2.17 


8.38 


7.16 


40 + long 




9,10 


LW20 


00 24 52.1 


-71 56 11 






8.31 


7.11 


49 




2,8 


LW21 


00 25 23.2 


-72 1 1 05 






8.44 


7.28 


38 




2,8 


LW22 


00 25 30.1 


-72 04 32 


12.06 


2.48 


8.16 


6.93 


63 + long 







Notes: Column 2 and 3 give coordinates from 2MASS. Column 4 is the mean V brightness. J and K colours are from 2MASS converted to the 
AAO system, if not stated otherwise. If there is an indication for a long secondary period exceeding the time span of the monitoring, the star is 
marked with '+long'. Uncertain values are marked with a colon. 

Remarks: (1) see Lebzelter et al. <2005> : (2) no K or V - 1^, not in the CMD (Fig.0; (3) currently only the longer period is visible; (4) 
photometry in the CMD from Fox )1982t - see text; (5) cuiTent periodicity unclear, see Lebzelter et al. <2005t for details; (6) close neighbour 
on 2MASS images, NIR fluxes and position possibly wrong; (7) not in 2MASS point source catalogue, coordinates calculated from our images, 
K magnitude from Origlia et al. <2002t ; (8) only on Mount Stromlo frames; (9) no Mount Stromlo data; (10) stars with near infrared excess 
detected by Origha et al. iMill ; (11) see text in Section|3]for comments on the periods. 



provide only near infrared finding charts of low spatial resolu- 
tion. Three stars from their list can be identified unequivocally 
with optical counterparts (stars number 1, 3 and 5 from Origlia 
et al.). All three are variables, namely LWIO, LW14 and LW19, 
respectively. The source number 4 from Origlia et al. is prob- 
ably related to the variable star LW7. The fifth star listed in 
Origlia et al. (their number 2) could not be identified with any 



LPV. It is the star with the weakest infrared excess in their list. 
It is surprising that the star (LWl 9) with the largest mass loss 
rate in Origlia's list (number 5) is neither one of the brightest 
variables nor has it a long period. A similar comment can be 
made about VI 8, the variable with the largest infrared excess 
in a survey of 47 Tuc by Ramdani & Jorissen J200II see also 
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0.0 0.2 0.4 0,6 0.S 



2.0 0,0 0.2 0,4 0.6 O.S 1.0 



2.0 0.0 0,2 0.4 0.6 O.K 1.0 




0.0 0.2 0.4 0.6 Q.S 1.0 



2.0 0.0 0.2 0.4 0.6 O.S 1.0 



2.0 0.0 02 0.4 0.6 O.S 1.0 



Fig. 2. Light curves versus phases for selected, monoperiodic variables from Table[T] The period given in the table was used to 
calculate the phase. 



Lebzelter et al. l2005l l. There is obviously no strict correlation 
between pulsation period and mass loss rate. 



5. The K-\ogP diagram 



Given the large number of LPVs now known in 47 Tuc, we are 
now in a position to examine the PL relations in the cluster 
We do this in the form of the K-log P relation, which is shown 
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V19 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



Fig. 3. Colour-magnitude diagram of 47 Tuc. Variable stars are 
indicated by open circles. Asterisks mark stars taken from Fox 
J1982> while the triangle is V3 using data from Eggen (I1975> . 
For all variables mean magnitudes and colours were used. 

in Figure|3 Generally, we do not have K light curves for the 
LPVs in 47 Tuc so we can not readily calculate mean K mag- 
nitudes. However, a small number of the variables have been 
observed many times (see Fox 119821 Menzies & Whitelock 
1 19851 and Frogel, Persson & Cohen [T^81> and good estimates 
of the mean K magnitude and K amplitude can be obtained. For 
each star plotted on Figure|4] the K magnitude is the mean of 
the maximum and minimum observed K magnitude. For stars 
without multiple observations, the 2MASS K magnitude has 
been used. All / and K magnitudes have been converted to the 
AAO system using the conversions in Allen & Cragg l20QjJ 
and Carpenter J200H : these are the values shown in Table [2 
A comparison of the K amplitude with the V light amplitude 
obtained in the present study shows that the K amplitude is ap- 
proximately 20% of the V amplitude. The error bars in Figure|4l 
have a full length of 20% of the full visual amplitude of the 
pulsation mode associated with each point. Hence, they should 
represent the K amplitude of the variables. 

Also shown in Figure|3 are the sequences of LPVs in the 
LMC as given by Ita et al. (.2_004j. The magnitudes have been 
adjusted to bring them from an LMC distance modulus of 18.55 
to a 47 Tuc distance modulus of 13.50. The large-amplitude 
Mira variables fall close to sequence C, while the bulk of the 
smaller amplitude variables fall near, but not necessarily on, se- 
quences B^, and C. The lack of variables on sequences A^ 
and A" may be real or it may be that these stars have amplitudes 
that are below our detection limit of ~0. 1 mag. 

It is notable in Figurel^that the 47 Tuc LPVs do not ap- 
pear to fall exactly on the LMC sequences. For example, the 
majority of the smaller amplitude variables brighter than K-1 
fall between sequences B^ and C while fainter than K-1, most 
variables fall between sequences A" and B" with a small num- 
ber appearing to fall on sequence C. This is probably a result 
of the different masses of the stars in the 47 Tuc and LMC sam- 
ples. We address this possibility in the next section. 
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Fig. 4. The K-\og P diagram for the LPVs in 47 Tuc. The error 
bars are scaled V light amplitudes. Only variables with rea- 
sonably well-determined periods (periods without a following 
colon in Table[0 are shown. The lines show the sequences of 
LPVs in the LMC as given by Ita et al. j2004l l. Sequence C is 
thought to be linked to fundamental mode pulsation, the other 
sequences to first to third overtone pulsation. 

If we neglect V19 which appears anomalously blue (see 
Sect. 13, stars seem to evolve up to a: 6.7 in the 47 Tuc equiv- 
alents of sequences B and C, and to then switch to sequence 
C for the final stage of evolution. Using the results in Wood 
et al. JI999> . this switch seems to correspond to a transition in 
pulsation from a low order radial overtone mode to the funda- 
mental radial mode (compare also Lebzelter et al. l2005t . 

The tip of the RGB in the K-\og P diagram occurs at K -1 . 1 
(see Fig.|6} at log L/L© =3.355. From Figs.0]and|5]we can say 
that about half of the variables are more luminous than the RGB 
tip and therefore must be AGB stars. The stars below the RGB 
tip could be either on the AGB or the RGB. 

6. Pulsation models 

6.1. Description of the models 

With a substantial set of LPVs now known in 47 Tuc, we 
are in a position to make theoretical models for the period- 
luminosity laws. For our models, we adopt the following pa- 
rameters for 47 Tuc (see the Introduction): distance modulus 
(ot - M)v-13.5, reddening E{B - y)=0.024, helium mass frac- 
tion Y=0.27, metal abundance Z=0.004 and main-sequence 
turnoff mass 0.9 Mq. Since we are computing pulsation peri- 
ods, it is important to get the radii (and hence Teff) of the mod- 
els correct. In order to realize this condition, the mixing-length 
in the convection theory was set so that the models coincided 
with the giant branch of 47 Tuc. Figure [S] shows the models in 



8 



T. Lebzelter and P.R. Wood: Long period variables in 47 Tuc: direct evidence for lost mass 



3.5 



2.5 





1 1 1 1 1 1 1 1 
■ ■ 


1 1 1 1 1 1 1 

■ 

■ □ 


- 






- 
- 












- 


■ • o 












o 








— Fiducial giant branch 


- 


° 


Variable stars 




Models: 




o 


RGB, no mass loss 




□ 


AGB, no mass loss 




• 


RGB, with mass loss 




■ 


AGB, with nnass loss 




, , 1 , , , 1 . 


, 1 , , , 1 , , , 



I I I I \ I I I I I I I I I I I I I I I I 

3.62 3.6 3.58 3.56 3.54 3.52 



Fig. 5. The tip of the 47 Tuc giant branch in the HR-diagram. K 
and J-K of the variables were converted to log L/Lq and log Teff 
using the transforms in Houdashelt et al. (12000 alEOOObi) . The 
fiducial giant branch given by Hesser et al. (I1987> is also 
shown, with V and B-V also converted to log L/Lq and log Teft 
using the transforms in Houdashelt et al. The pulsation models 
are plotted using symbols identified on the figure. 



the HR-diagram. The giant branch of the models clearly passes 
through the region of the giant branch occupied by the vari- 
ables. 

The properties of the models are given in Table |2] The lin- 
ear non-adiabatic pulsation models were created with the pul- 
sation code described in Fox & Wood tl982i . updated to in- 
clude interior opacities of Iglesias & Rogers J1993> and low 
temperature opacities of Alexander & Ferguson ( I1994> . This 
code uses a mixing length theory of convection that explicitly 
treats variation of the convective velocity with time. The core 
mass Mc was obtained from the L-M^ relation of Boothroyd & 
Sackmann (,1988 | . The most luminous RGB model corresponds 
to the RGB tip luminosity of the 0.9 Mq, Z=0.004 tracks of 
Fagotto et al. ( I1994> . 

The models with mass loss were constructed according to 
the following prescriptions. Firstly, note that without mass loss 
the giant branch should extend to very high luminosities well 
beyond logL/L© = lO'^L/L©. Since the most luminous stars on 
the 47 Tuc giant branch have L/L© x 4OOOL/L0, it is assumed 
that mass loss dissipates the hydrogen-rich envelope at about 
this luminosity. Evolutionary models with mass loss were con- 
structed adopting a Reimers' mass loss law (Reimers I1975> . 
and the evolution rates on the the giant branch given by the 
evolutionary tracks of Fagotto et al. ( 119941 . The Reimers mass 
loss rate was multiplied by a factor j] which was adjusted so 
that the stars left the AGB at L/Lq 4000 L/Lq (a value 77 = 
0.33 was required). The resulting models are shown in the HR- 



diagram Figure |5] and the masses of the models are given in 
Table|2l 

The K-logP diagram for models both with and without 
mass loss is compared with the observed K-logP sequences 
in Figure |6l We now consider the small and large amplitude 
sequences separately. 

6.2. Small amplitude variables 

It is clear that the models without mass loss fail to reproduce 
the observed periods for the smaller amplitude pulsators (the 
overtone pulsators, which have log P < 2). In fact the model 
sequences avoid the observed sequences. Since the model cal- 
culations are for linear pulsation models, they should fit the 
small ampUtude variable sequences. 

A simple way to get the observed and model sequences to 
agree would be to make the giant branch cooler at a given lu- 
minosity. Since the pulsation period for the overtones P cc Ri 
(Fox & Wood [T982> . and since a given lu- 

minosity P oc r^jf3. In order to fit the theoretical overtone se- 
quences to the observed sequences, log P needs to increase by 
~0.15, which means that logTeif needs to decrease by ~0.05. 
As can be seen from Figure|5l this is far greater than any likely 
uncertainty in log T^ff. We therefore conclude that models that 
have not undergone considerable mass loss since the main- 
sequence can not explain the pulsation periods of the overtone 
LPVs in 47 Tuc. 

In contrast to the models without mass loss, the overtone 
pulsation periods for the stars that have undergone mass loss 
agree well with the observations. The periods of the mass loss 
models are longer than those of the models without mass loss 
due to the lower stellar mass (P oc M 2), and the slope of the 
A'-log P relations is smaller due to the decrease in mass with 
luminosity. The decreased slope is particularly prominent for 
the second overtone which moves from having a period close 
to that of the 3rd overtone at low luminosities to a period close 
to that of the 1st overtone at high luminosities. 

6.3. Large amplitude variables 

In contrast to the small amplitude variables, the large- 
amplitude variables i.e. the Miras are consistent with theoreti- 
cal models without mass loss, while the models with mass loss 
do not fit the observed K-logP sequence. (We ignore V19 in 
this comparison since Figure |3l suggests this star is a Galactic 
halo star rather than a 47 Tuc star.) 

We believe this contradiction can be explained by the non- 
linear effects in the pulsation of red giants with large ampli- 
tudes. Some preliminary nonlinear pulsation calculations made 
for these models show that the full-amplitude pulsation peri- 
ods are considerably shorter than the linear periods due to a 
change in the envelope structure associated with large ampli- 
tude pulsation (see also Ya' Ari & Tuchman I 19961 1. An example 
of the change in period of one of these models, perturbed with 
a low amplitude and allowed to reach limiting amplitude, is 
shown in Figure^ Thus the periods of nonlinear fundamental 
mode models may be able to explain the periods of the Miras in 



T. Lebzelter and P.R. Wood: Long period variables in 47 Tuc: direct evidence for lost mass 9 



Table 2. The pulsation models 
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3.5949 
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0.8166 
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1.90 
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49.5 


30.5 
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17.3 
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0.4417 
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69.2 
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3.5629 
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1.90 
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53.2 
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Notes: {[Hp is the ratio of mixing-length to pressure scale height. Pq, Pi , P2 and P3 
are the linear periods (in days) of the fundamental, 1st, 2nd and 3rd overtone modes. 



47 Tuc, although this needs further exploration. The nonlinear 
models will be described elsewhere. 

6.4. Evolution of pulsation mode with luminosity 

Figure |6| provides an indication of the pulsation evolution of 
stars evolving up the giant branch. Ignoring V19 (as noted 
above), it seems that stars start pulsating at K ~ 7.7 in the 
1st to 3rd overtone, then evolve up to higher luminosities to 
K ~ 6.9 where the stars transit to fundamental mode pulsation. 
Further evolution is in the fundamental mode, until mass loss 
terminates AGB evolution at K ~ 6.2. This is broadly consis- 
tent with what is expected from the linear non-adiabatic growth 
rates of the models. These are shown plotted against luminosity 
in FigurelH It should be remembered that the theoretical growth 
rates for these highly convective stars are very uncertain and 
should be regarded as indicative only. The general behaviour 
shown in Figure|8lis that, at low luminosities the 2nd overtone 



has the highest growth rate, at intermediate luminosities the 1st 
overtone has the highest growth rate, and at high luminosities 
the fundamental mode grows most rapidly. The 2nd overtone 
also has a high growth rate at high luminosities but it is likely 
to be overwhelmed by the fundamental mode in the nonlinear 
case. This suggests an evolution from 2nd to 1st overtone and 
then to fundamental mode. This is similar to the observed sit- 
uation, although the overtones co-exist rather than following a 
distinct progression with luminosity. 

7. Summary and Conclusions 

We have identified 22 new variable red giants in 47 Tuc and de- 
termined periods for another 8 previously known variables. All 
red giants redder than V-Ic = 1.8 are variable at the limits of our 
detection threshold, which corresponds to 6V~Q.l mag. This 
colour limit corresponds to a luminosity logL/L0=3.15 and it 
is considerably below the tip of the RGB at log L/L0=3.35. In 
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log P(days) 

Fig. 6. The K-log P diagram for models both with and without 
mass loss. The fundamental mode and the first three overtones 
are shown. The variables with reasonably well-determined pe- 
riods are also shown, as in Figure^] 



the K-logP diagram, the 47 Tuc variables do not closely fol- 
low the ridge lines of PL relations seen for LPVs in the MCs, 
indicating that the PL relations are mass dependent. 

Linear non-adiabatic modelling was used to try to repro- 
duce the observed PL relations, especially for the low ampli- 
tude pulsators where linear calculations should be appropri- 
ate. It was shown that models without mass loss can not re- 
produce the observed PL laws for the low amplitude pulsators. 
Models that lose sufficient mass to terminate AGB evolution 
near L ~ 4OOOL/L0 do reproduce the observed PL relations 
for low amplitude variables. This is the first time that measure- 
ments of the masses of stars on the AGB have shown that mass 
loss of the order of 0.3 occurs along the RGB and AGB. 

The linear pulsation periods do not agree well with the ob- 
served periods of the large amplitude Mira variables, which 
pulsate in the fundamental mode. The solution to this prob- 
lem appears to be that the nonlinear pulsation periods in these 
low mass stars are considerably shorter than the linear pulsa- 
tion periods due to a rearrangement of stellar structure caused 
by the pulsation. Although such effects have been seen in pul- 
sation models before, the 47 Tuc stars studied here provide the 
first observational evidence for this effect. 

The observations show that stars evolve up the RGB and 
first part of the AGB pulsating in low order overtone modes. 
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Fig. 7. Luminosity L, pulsation period P and radii of several 
mass points plotted against time in a model with M - 0.6567 
Mo and L = 3 162 L©. The linear pulsation period of this model 
is 240 days while the nonlinear pulsation period after 67 years 
of relaxation is 185 days. 



a" 0.004 



-0.002 




3.2 3.4 

log L/Lg 



Fig. 8. The growth rate to^ (days"' plotted against log L/L© for 
the models in the mass loss sequence. Wi is the real part of 
the eigenvalue <y where a time dependence exp(wf) is assumed. 
The K magnitude for AGB stars is shown on the top axis. 

then switch to fundamental mode at high luminosities. The lin- 
ear non-adiabatic growth rates of models suggest that such be- 
haviour should occur but the models at this stage are only in- 
dicative. It is hoped that future improved models including the 
effect of turbulent viscosity (e.g. Olivier & Wood 2005 1 will 
allow a reliable determination of the growth rates and mode se- 
lection processes in red giant stars, as well as an estimation of 
the eff'ect of nonlinear pulsation on the pulsation period. 
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